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Abstract 



The gravitational-radiation-induced inspiral of a binary system of compact objects 
is considered. A scheme is described to model the regime in which the gravitational 
interaction is too strong to use weak-field approximation methods, but the time scale 
for decay of the orbits is still long compared to the orbital period, by numerically solving 
for a stationary spacetime which approximates the slowly evolving one. Equilibrium 
is to be maintained in the radiating system by imposing a balance of incoming and 
outgoing radiation at large distances. Numerical results from non-linear scalar field 
theory have shown that such an approach can be effective modelling a slowly evolving 
solution to a wave equation. 
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1 Introduction 

A pair of compact objects (black holes or neutron stars) in binary orbit will, according to 
Newtonian gravity, remain in the same orbits forever. In general relativity, however, the 
system will emit gravitational radiation, causing the bodies to spiral in towards one another. 
The gravitational radiation given off by this system is a prime candidate for detection by 
upcoming gravitational wave telescopes such as VIRGO and LIGO [|TJ. 

Full three-plus-one numerical evolution of the Einstein equations is a powerful tool for 
determining the history of such a system and the form of the gravitational radiation emitted, 
but is limited in its applicability by instabilities which prevent the simulation from being 
evolved for times longer than several orbital periods. Thus it is not only efficient but also 
necessary to limit full-numerical simulations to the portion of the evolution for which there is 
not an applicable approximation scheme [fj|. The early stages of inspiral can be handled by 
post-Newtonian methods |§ , while the final post-merger ringdown can be treated with black 
hole perturbation theory Q . The purpose of the quasi-stationary approximation is to provide 
improved early waveforms and later-time supercomputer initial conditions by modelling a 
late inspiral phase, in which some non-perturbative gravitational effects are relevant, but the 
radiation-induced inspiral is still occurring slowly. 

2 General Formulation 

2.1 Quasi-Stationary Approximation 

The idea, initially proposed by Steven Detweiler || ||, is that if the inspiral is slow, the system 
is nearly periodic: after one orbit, the objects have returned almost to their original locations, 
and radiation which has moved out has been replaced with new radiation of approximately 
the same shape. If the objects' orbits are circular rather than elliptical, the spacetime is 
nearly stationary. If the approximate orbital frequency is Q, moving forward in time by 5t 
and rotating the resulting spatial slice by —Q 5t will not change the picture very much. 

This approximation can be used to simplify the numerical problem by solving for an 
exactly stationary spacetime which serves, over some period of time, as a reasonable approx- 
imation for the slowly evolving spacetime. In the process, the three-plus-one dimensional 
evolution problem is reduced to a three-dimensional, non-dynamical one, not only reducing 
the size of the computational grid, but also hopefully avoiding evolutionary instabilities. 

Since gravitational radiation emitted from the orbiting system carries away energy, some 
modification must be made to the physical problem in order to allow this equilibrium solution. 
One approach fA, RL M is to require the spatial geometry of the stationary spacetime to be 
conformally flat and solve for the conformal factor using only a subset (the constraints) of 
the Einstein equations. This method enforces stationarity by in effect discarding degrees of 
freedom associated with the radiation which would carry away energy. Our approach, on 
the other hand, is to keep the gravitational radiation, but nullify the radiation reaction by 
balancing the outgoing radiation by an equal amount of incoming radiation. In so doing, 
we will solve the full Einstein equations in the presence of the actual sources, and simply 
replace the physical boundary condition of outgoing radiation at large distances with one of 
balanced radiation. 
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2.2 Radiation- Balanced Boundary Conditions 

In the general context of a theory in which a pair of orbiting sources for a wave equation 
give off radiation, it is instructive to consider three types of solutions. In a Type I solu- 
tion, the radiation is outgoing, carrying away energy and causing the orbits to decay In 
a Type II solution, the radiation is again outgoing, but some additional force, uncoupled 
to the radiating field, keeps the sources in their circular orbits, resulting in a stationary 
spacetime. A Type III solution is also stationary, with the sources remaining in circular 
orbits of a constant frequency, but this time there is no external force, and the equilibrium is 
maintained by a balance of incoming and outgoing radiation. Type I is the physical situation 
we ultimately wish to model. Type II should be a reasonable approximation to Type I if 
the Type I solution in question is inspiralling only slowly. It should not be appropriate to 
general relativity, in which all matter and energy couples to the gravitational field, but is 
useful for comparison in theories which allow for external forces. Type III is the stationary 
solution we wish to find numerically, and, where appropriate, to relate it to a Type II or I 
solution. 

To analyze in detail the implications of boundary conditions on radiating, stationary so- 
lutions (Types II and III), we have considered the theory of a non-linear scalar field ip(t, p, 4>) 
in a 2 + 1 dimensional Minkowski spacetime, the simplest theory with nonlinearity, orbits, 
and radiation. If the source is required to rotate at a constant angular velocity Q, so that 
the charge density a(t, p, (ft) is a function only of p and tp = <fi — Qt, there exist stationary 
solutions for which the field exhibits the same symmetry; in that case, the wave equation 
can be written 



where p) is a non-linear coupling. A convenient choice of charge distribution a is a 
positive charge at a radius a and azimuthal angle p = and an equal and opposite charge 
at <p — 7r and the same radius. This distribution satisfies cr(p, — <p) = cr(p, (p). 

In the context of a numerical solution on a finite grid, the Type II (purely outgoing 
radiation) solution ip^ ut is defined by applying a Sommerfeld (ip^+c^ipj = V'.p - ^ c — 0) 
boundary condition at a the radius R of the outer boundary of the grid. If the non-linear 
term becomes small at large distances, the form of any solution near a large- R boundary will 
be a solution to the vacuum, linear version of ([!]); each angular Fourier mode will be a linear 
combination of two independent solutions, which can be identified as a purely ingoing and 
a purely outgoing solution. The limit ip out : = lim^oo ip^ nt is simply the solution consisting 
only of outgoing modes, with the coefficients of all the ingoing modes set to zero. [The 
entire discussion can be repeated for ingoing radiation, with the resulting solution being 

The most familiar Type III radiation-balanced (RB) solutions are standing-wave solutions 
in which the Dirichlet (ip = 0) or Neumann (ip tP = 0) boundary condition is enforced at p = 
R. If the large-distance form of any of these solutions is resolved in angular Fourier modes, 
it is found that the amplitudes of the coefficients of the ingoing and outgoing components 
are equal; however the relative phase of the two independent vacuum solutions depends on 
the choice of R, as does the radiation amplitude corresponding to a given source strength. 
Unlike outgoing- (or ingoing-) wave solutions, standing-wave solutions do not tend towards 
any limit as the boundary radius is taken to infinity. 
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Standing waves are not the only solutions with equal magnitudes of incoming and outgo- 
ing radiation; in the case of a linear theory (A = 0) another such solution is simply a linear 
superposition of the ingoing and outgoing solutions (LSIO) t/'lsio = | (^in+Vw)- For a given 
source strength, this turns out to be the RB solution with the smallest amplitude, which 
means in effect that it's the solution with just the radiation needed to keep the sources in 
equilibrium and no "extra" radiation. From the point of view of using a Type III solution to 
approximate a Type II one, we can take this "minimum energy radiation balance" (MERB) 
solution, extract its outgoing component, and identify it with the outgoing radiation. 

In the non-linear theory, the superposition of two solutions is no longer a solution, but 
the same effect can be achieved by reformulating (HD as a Green's function problem and using 
the average of the advanced and retarded Green's functions to find a MERB solution. We 
have done this |L(J and found that even for highly non-linear theories, this time-symmetric 



Green's function method can be used to find a good approximation for the outgoing solution. 
This lends confidence that in the case of general relativity, if we can find a Type III solution 
with little or no "superfluous" radiation, we may be able to relate it to the physical Type I 
solution in an analogous way. 

As an aside, it is worth stressing that while the equation (|]) is elliptical inside and hy- 
perbolic outside the "speed of light circle" p = 1/fi, we experienced no difficulties finding 
the numerical solution to the problem using closed-surface boundary conditions usually as- 
sociated with a purely elliptical equation. In particular there were no discernable artefacts 
at the light circle, which required no special treatment in the numerical solution. 



3 Cosmic String Toy Model 

To apply the quasi-stationary method to general relativity, we first note that the stationarity 
of the numerical spacetime is described by a Killing vector K ~ d t + Qd^. This means 
that the geometry of the four- dimensional spacetime is determined by the geometry on a 
representative three-dimensional surface. Before attacking this three-dimensional problem, 
we first consider a computationally simpler toy model in which the sources and fields posess 
an additional Killing symmetry K\ ~ d z corresponding to a translation perpendicular to 
the orbital plane. This then is a model describing not orbiting localized sources, but rather 
orbiting infinitely long lines of mass, or cosmic strings.^] The addition of a second Killing 
vector field in four- dimensional spacetime means that the geometry need only be numerically 
determined on a two-dimensional surface. Our primary intention in looking at this model is 
to test the method rather than to obtain astrophysically meaningful results, but it's worth 
pointing out that this would provide a solution to the two-body problem in full vacuum 
general relativity, albeit with somewhat unusual boundary conditions. 

While numerically simpler, the problem of orbiting cosmic strings brings with it some 
conceptual subtleties, which are described in the next two subsections. 



lr The cosmic strings in question posess a curvature singularity along the string and generate Riemann 
curvature in the external spacetime. They are not to be confused with the more familiar special case of 
cosmic strings consisting of a conical singularity surrounded by a flat external spacetime. 
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3.1 Spacetimes with Two Killing Vectors 

One question is how to choose a basis which takes full advantage of the two Killing symmetries 
of the stationary orbiting cosmic string spacetime, thereby removing any gauge freedom from 
the numerical problem. Since the orbits of K (which is associated with the stationarity of 
the spacetime and also present in the three-dimensional case) are not surface forming, a 
coordinate basis on the spacetime four-manifold turns out not to be ideally suited to the 
task. Instead, one defines coordinates x 2 and x 3 on the two-manifold of Killing-vector orbits 
(every point of which represents an equivalence class of points in the four-manifold connected 
by Killing trajectories. The associated basis vectors e<i and e% commute when restricted to 
the two-manifold, but the corresponding vectors they define on the four-manifold do not 
commute. Thus ei and e%, together with the Killing vectors K and Ki, define a non- 
coordinate basis on the spacetime four-manifold. 



Extending the results of of [fL^l , we have used this method to find explicit forms of the 



Einstein equations, both in a general two-Killing-vector spacetime and in the case of orbiting 



cosmic strings |TT| . The most convenient choice of coordinates (p, <p) on the two- manifold of 



Killing-vector orbits are the anologue of the corotating polar coordinates used in Section [L2 



3.2 Perturbations to Levi-Civita 

To analyze and apply boundary conditions to the radiation far from the sources, we need to 
cast the spacetime in a form consisting of a static background plus a perturbation describ- 
ing gravitational radiation on that background. With localized sources, that background 
could in principle be taken to be Minkowski spacetime, but this is impossible with line-like 
sources because the spacetime of an infinite line of mass is not asymptotically flat, even 
as one approaches infinity in directions orthogonal to the line. (Equivalent ways of seeing 
this are to note that the gravitational potential of a line mass falls off logarithmically, or 
that the quantity mass-per-unit-length is dimensionless in general relativity.) The simplest 
appropriate background spacetime is the static, cylindrically symmetric Levi-Civita space- 
time [|13| describing a single cosmic string; the Riemann curvature is determined by a single 
dimensionless parameter C which is associated with the mass per unit length of the string. 



Analysis of perturbations to Levi-Civita shows that tensor as well as scalar modes 
are governed by an equation 

& p i> m (p) + P~ l d p ijj m {p) + F m {p)tfj m {p) = 0. (2) 

The potential F m (p) depends on the tensor character and parity of the mode considered, but 
the leading term at large p is 

F(p) ~ m 2 n 2 p- 2C ^- c+c2 \ (3) 

[Note that for the flat background case C = 0, this reduces to the large- p limit of the linear 
version of (0).] 

3.3 Strategy for Radiation Balance in General Relativity 

The general solution to (0) is, as before, a linear combination of an outgoing and an ingoing 
solution. In principle, a general RB solution can be defined as one with equal ingoing and 
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outgoing amplitudes. However, for the purpose of numerically enforcing RB, it is simpler 
to use an equivalent condition. If the locations of the sources are taken to be at ip = and 
if = 7r, the condition of RB is equivalent to invariance of the geometry under the reflection 
if — > — f. This can be seen by noting that the flux T tp of energy through a surface of constant 
p will be proportional to i/) tP ip tip . If ip is even under <p- inversion, ip jP will be even and ip jV> odd, 
making T tp odd under f — > —f>. This means that the total rate of energy (integrated over <fr 
or equivalently ip) flowing through a constant-p surface will vanish for a field configuration 
obeying this reflection symmetry. If the sources are taken to have equal mass-per-unit-length, 
the geometry will also be unchanged by a 180° rotation which interchanges them, meaning 
that one need only solve for the metric over the interval ip G [0, 7r/2]. This "quadrant 
symmetry" is the analogue in the two-dimensional problem of the octant symmetry familiar 
from |7], [5], fj; it is not only a computational simplification but also a means of automatically 
ensuring a balance of radiation. 

As noted in Section [2.2| , a wave equation has not one but an entire family of RB (Type III) 



solutions. We would like to choose the one which can be most directly related to the (non- 
stationary) Type I solution with outgoing radiation. While the nature of the full nonlinear 
Einstein equations means that the Green's function method used for scalar field theory will 
not work in the gravitational case, we can still try to select a RB solution with as little 
gravitational radiation as necessary, applying the following scheme: 

1. Treating the sources perturbatively, obtain a guess for the relative phase of the ingoing 
and outgoing components of the radiation and thus a candidate asymptotic form for 
each mode of the MERB; 

2. Apply a boundary condition to the values of the metric components at the outer- 
most two p values in the grid which forces the perturbation to be proportional to the 
candidate MERB; 

3. Modify the candidate MERB and repeat step 2 to get a family of solutions, and find 
the one with the smallest total gravitational wave amplititude. 



4 Conclusions and Outlook 

This paper has described an ongoing research program to find numerically a stationary 
spacetime which can approximate over some stretch of time a slowly inspiralling compact 
object binary. A stationary solution to the radiative problem is to be achieved by replacing 
the physical boundary condition of purely outgoing radiation with a condition corresponding 
to an equal mix of ingoing and outgoing radiation. 

Thus far, the numerical work has been limited to testing the usefulness of this sort 
of radiation-balanced approach in non-linear scalar field theory, but we have formulated a 
toy problem in general relativity with infinitely extended sources and are about to begin 
analyzing it numerically. Looking beyond, to the application of the method to the physical 
problem of localized sources, it is instructive to consider how such an equilibrium three- 
dimensional spacetime could be used. 

First, since the spacetime will contain at large distances a superposition of ingoing and 
outgoing gravitational radiation, we should be able to separate out the outgoing-wave con- 
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tribution and use it as an approximation for the purely-outgoing radiation from the physical 
system. 

Second, a section of the spacetime could be used as an alternative to post-Newtonian or 
conformally flat initial data for a full numerical evolution. (Since our approach contains in 
some sense "too much" and the conformally flat approach "too little" radiation, a comparison 
of the two should provide insight into the effects of the radiation on the evolution of the 
spacetime.) It should be pointed out that while the RB spacetime is an equilibrium solution 
of Einstein's equations, valid for all time, it can still be used to provide initial data for a full 
numerical evolution, since the outer boundary condition of such an evolution will be one of 
outgoing rather than balanced radiation, and thus once the "missing" ingoing radiation has 
failed to propagate from the outer boundary to the location of the sources, the orbits will 
begin to decay. (The author is indebted to B. Briigmann for clarifying this point.) 
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